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Abstract 

A heuristic but pedagogical derivation is given of an explicit for- 
mula which accurately reproduces the period of a simple pendulum 
even for large amplitudes. The formula is compared with others in 
the literature. 



1 Introduction 

The dynamics of a simple pendulum for small amplitudes is probably the 
most widely used example by educators to illustrate simple harmonic motion: 
Both the experiment and theory are accesible to college students. However 
the discussion of the large amplitude oscillations is rarely carried out, prob- 
ably because there is no "simple" analytical formula for the period. Over 
the years several approximation schemes have been developed to discuss the 
large angle oscillations and in this paper we provide yet another. It turns 
out that the formula derived here, which is quite accurate even for very large 
amplitudes, was obtained a long time ago by Ganley jl] using a very differ- 
ent method. However, while the final result is the same, the derivation here 
is probably simpler and also directly linked to the usual derivation for the 
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linearised pendulum. Indeed, the interested reader could potentially use the 
presentation here as the starting point for a systematic improvement of the 
formula. Furthermore, an important difference between the result of Ganley 
and that obtained in this paper is mentioned at the end of the next section. 



2 Approximation Scheme 

The equation of motion for the simple pendulum can be derived by applying 
Newton's Second Law to motion along the tangential direction of the arc, 
giving 

+ ysin0 = O, (1) 

where / is the length of the pendulum, g the acceleration due to gravity and 
9 the angular displacemnent. The usual linearised approximation involves 
replacing sin^ by 9. Let g/l — cUq, so that the period, T , of a linear pendulum 
is T = 27i-(c<Jq)~ 1 / 2 = ^, which up to a constant of proportionality is the 
unique expression obtained by dimensional analysis. Indeed, by dimensional 
analysis the period of the non-linear pendulum must be of the form 

T a = T f(a) , (2) 

where a is the amplitude of oscillations and f(a) a function to be determined. 
An approximate expression for f(a) can be obtained by first rewriting the 
equation of motion in the suggestive form 

s+ ^(^fiy =0 . (3) 

If one could treat the term in brackets as a constant, as it approximately 
is for small angles, then an attempted estimate of the period of the non-linear 

pendulum would be T a « 2vr (co^-^f 1 ) . Of course this formula does not 
make much sense as it stands since the period cannot depend on the angular 
displacement: The varying 9 should be replaced by some "average" value f3 
that depends on the amplitude a. Thus we obtain 

t . . /sin(/?)\~ 1/2 

/(«) ~ [-p) ■ ( 4 ) 
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Perturbative solutions [2] of the non-linear pendulum equation of motion 
give f{u) ~ 1 + Yg. Comparing this with the perturbative expansion of the 
right-hand-side of Eq. (@J implies 

16 * 12 ' (5j 

or (3 w ^l 2 . Hence an first approximation for the period of a non-linear 
pendulum is 

The quality of this approximation can be determined by comparing it to the 
exact value given by 

2T Q T ^ ( . 2 a, 



T ex = ~^K ^in 2 (^)j , (7) 
where 

K(m) = f /2 , (8) 

jo [1 — m sin y) L ' 1 

is the complete elliptic integral of the first kind. The two expressions T al /T 
and T ex /T are plotted in Figure(l) which shows that the approximation T al 
gives an accuracy of 1% for amplitudes as large as a ~ 2.25 radians. It 
should be noted that the approximation always lies below the exact answer, 
but unlike the exact answer it does not diverge at a = it. 

While the same formula for the ratio T al /T was obtained by Ganley in 
an interesting paper pQ using different methods, the approach presented here 
gives T a \ explicitly and consequently in the limit of small angles the known 
expression T is also recovered. 



3 Comparison with Other Approximations 

Perturbative solutions for f(a) are known (or can be obtained by expansing 
the right-hand-side of (jZJ)). From the symmetry of the problem, one sees that 
only even powers of a can appear. A natural question is: "How many terms 
of the series expansion must be kept to get a 1% accuracy in the period for 
an amplitude a ~ 2.25?" (that is, to compare with the approximate formula 
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of the last section). A simple numerical exercise shows that all terms up to 
and including the eighth-order term must be taken into account. 

In a recent paper Kidd and Fogg [3j motivated the approximation f(a) ~ 



l/ycos(a/2) which provides a 1% accuracy in the period for a up to ir/2. 
Though the expression of Kidd and Fogg is accurate over a smaller range of 
amplitudes than the expression T al , it shares one interesting feature with the 
exact expression: A divergence at a = n. This suggests that one generalise 
the method of the last section by trying an ansatzior f(a) that also diverges 
at the extreme amplitude a = ir. One simple possibility, suggested by the 
form of Eq. (j3J) is 



where 7 is a positive constant to be determined. By comparing the first 
nontrivial terms of the power-series of both sides of the last equation, one 
obtains 7 = — 3/8. The resulting expression 



always exceeds the exact answer and has an accuracy of 1% for amplitudes 
up to a ~ 1.9 radians. Interestingly, the result T a2 is precisely the expression 
of Molina jl] but obtained there with a slightly different motivation. 

4 Conclusion 

The approach presented in Sect. (2) may possibly be viewed as a kind of "adi- 
abatic approximation" and perhaps improved on to give successively better 
approximations. The main result (jHJ) is relatively simple and accurate to very 
large amplitudes. While the ratio ^ was previously obtained in pQ using 
different arguments, here an explicit expression for T a \ has been obtained. 
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Figure Caption 

Figure 1: Plot of T ex /T and T al /T against the angular amplitude a. 
The exact expression of course diverges for a = ir. 




6 



